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We present a method for studying the excitations of low-dimensional quantum spin systems based 
on the Jordan- Wigner transformation. Using an extended RPA-scheme we calculate the correlation 
function of neighboring spin flips which well approximates the optical conductivity of Sr2CuC>3. We 
extend this approach to the two-leg S = —ladder by numbering the spin operators in a meander-like 
sequence. We obtain good agreement with the optical conductivity of the spin ladder compound 
(La,Ca)i4Cu2404i for polarization along the rungs. For polarization along the legs higher order cor- 
relations are important to explain the weight of high-energy continuum excitations and we estimate 
the contribution of 4- and 6-fermion processes. 

PACS numbers: 75.10.Jm, 75.40. Gb, 75.40. Mg, 74.72.Jt, 75.30.Et 



I. INTRODUCTION 

Due to the presence of strong quantum fluctuations, 
low-dimensional spin systems show very complex behav- 
ior and provide a challenge for theoretical treatments. In 
this context, S= ^-Heisenberg spin ladders are especially 
interesting because they represent an interjacent system, 
in between the antiferromagnetic S= ^-Heisenberg chain 
and the two-dimensional antiferromagnetic Heisenberg 
model. Early on, these spin ladders were considered as 
systems which display a dimensional crossover between 
one and two dimensions. 1 However, spin ladders do not 
constitute a "smooth crossover", because even-leg spin 
ladders have a spin liquid ground state and finite spin 
gap. 2 ' 3 This is in contrast to the critical systems, the 
S=i-chain and the odd- leg ladders, which have alge- 
braically decaying spin correlations and to the 2D anti- 
ferromagnetic Heisenberg model, for which a long-range 
Neel-ordered ground state was established. 4 In a field the- 
oretical mapping of the low-energy modes on a 0(3) non- 
linear a model this even-odd effect in the number of legs 
results from the addition of a topological term. The term 
was demonstrated to be zero for even-leg ladders and the 
two-dimensional Heisenberg model and finite for odd-leg 
ladders as well as the S= ^ -chain. 5-8 Correspondingly the 
system is gapless for odd-leg ladders and gapful for even- 
leg ladders whereby the gap decreases exponentially with 
the number of legs. 

In this article we will focus on the antiferromagnetic 
two- leg S= ^-ladders. The two-leg ladder can be ap- 
proached conceptionally from the limit of strong coupling 
J_l along the rungs. 9-12 Then the elementary excitations 
can be considered as excitations of rung-triplets which 
propagate throughout the ladder due to the finite leg 
coupling J. For small coupling Jj_/J, a more natural 
description would seem to be in terms of the spinon ex- 
citations of the isolated legs. However, the excitations of 
the two-leg spin-ladders cannot be constructed perturba- 
tively from spinons of the chains since the rung coupling 
is a relevant perturbation. The spinons are confined and 
have to form bound states in the ladder-case. 13,14 



Of particular interest is the intermediate coupling 
regime J « J±, as this case is related to the two- 
dimensional systems. It is also realized in cuprate spin 
ladder compounds, 15 ~ 17 which are of interest due to their 
affinity to the cuprate high-T c -superconductors. This in- 
termediate regime is difficult to obtain and requires te- 
dious calculations. 

Here we propose an alternative approach, which is 
based on the Jordan- Wigner transformation 18,19 and 
which is fairly simple. In contrast to previous applica- 
tions of the Jordan- Wigner transformation to spin lad- 
ders by Dai and Su 20 and Azzouz et al. 21 we do not treat 
the phase factor within a mean-field or a flux-phase ap- 
proximation. Rather we expand the phase factor and 
treat the resulting interaction terms on the same foot- 
ing as the Ising term, which corresponds to a 4-fermion 
term in the Jordan- Wigner representation. We also apply 
this approach to the calculation of dynamic correlation 
functions. In particular we discuss the optical conductiv- 
ity cr(u>), as it provides valuable information about high- 
energy S'=0-spin excitations. In cr(u>) magnetic excita- 
tions can be observed due to the simultaneous excitation 
of a phonon, which breaks the inversion symmetry be- 
tween two neighboring spins. The mechanism, which has 
been suggested by Lorenzana and Sawatzky, 22 provides 
a finite dipole moment to a simultaneous flip of neigh- 
boring spins. Since the phonon can take arbitrary mo- 
mentum in this process, the magnetic S'=0-excitations 
have to be averaged over the whole Brillouin zone to 
ensure zero total momentum. This mechanism has re- 
cently been used to investigate the high-energy spin ex- 
citations in the spin-ladder compound (La,Ca)i4Cu2404i 
by Windt et al. 16 Using a perturbative approach, based 
on a continuous-unitary-transformation 23 and a Jordan- 
Wigner representation it was possible to identify the 
S'=0-bound state of two triplets in the optical conduc- 
tivity of (La,Ca)i4Cu2404i. 

Here, we discuss in more detail the Jordan- Wigner 
approach which we have previously used in Ref. 16. 
We introduce our approximation scheme for the Jordan- 
Wigner fermions for the simpler case of a Id spin-chain in 
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Sec. II. We calculate cr{u>) for the ld-spin chain in RPA 
and compare with the optical conductivity of Sr2CuO,3 
measured by Suzuura et al. 24 In Sec. Ill we generalize our 
Jordan- Wigncr treatment to two-leg ladders by arrang- 
ing the spins in a meander-like sequence. The mean field 
treatment of this representation, which is based on the 
expansion of all phase factors, is discussed in Sec. Ill A 
and the RPA-treatment in Sec. IIIB. We calculate dy- 
namic spin-flip correlation functions in Sec. Ill C and 
compare cr(u>) for the two-leg ladder with the optical con- 
ductivity of (La,Ca)i4Cu 2 404i in Sec. HID. The weight 
of processes involving the excitation of more than two 
fermions is estimated in Sec. Ill E. In Appendix A it is 
shown that our mean field treatment reproduces the cor- 
rect strong coupling limit, in Appendix B the full set of 
RPA-equations for the two-leg ladder is listed. In Ap- 
pendix C the role of the phase factor is explored, and 
in Appendix D the consequences of the artificial symme- 
try breaking introduced by our mean field treatment arc 
discussed. 



II. JORDAN- WIGNER TRANSFORMATION FOR 
THE ID SPIN CHAIN 

First we recall the Jordan- Wigner transformation for 
the Id spin chain and introduce our approximation 
scheme. The fact, that we find fairly good agreement 
with the spinon evaluation for the Id spin-chain, inspired 
us to extend the approach to two-leg spin-ladders in Sec- 
tion III. 

As spin operators do not obey canonical commutation 
relations, they are often transformed into either bosonic 
or fermionic operators, in order to permit the application 
of standard diagrammatic perturbation theory. With any 
mapping, however, the algebra of the original spin oper- 
ators has to be preserved. For the Jordan- Wigner trans- 
formation this is provided by rewriting the spin operators 
as fermionic operators with a long-ranged phase factor 

Sr =Ci e^, &=7r^ C t Cj , S? = (eje, - ±) . (1) 

j<i 

The representation is particularly useful for models 
with nearest-neighbor exchange interactions because in 
products of neighboring spins the phase factors drop out. 
This is due to the fact that the fermion operator c$ com- 
mutes with the phase of the same site ([cj, <pi] = 0) . 
Therefore the Id antifcrromagnctic Hciscnbcrg Hamilto- 
nian transforms to 

+ fa - \) (cl+iCi+i - } • (2) 

The first term in Eq. (2) corresponds to the XY-part 
of the original Heisenberg Hamiltonian. In the fermionic 



representation it acts as a nearest neighbor hopping. The 
second term, which originates from the Ising term of the 
original Heisenberg Hamiltonian, introduces a nearest 
neighbor density-density interaction among the fermions. 

A. Mean field treatment — d=l 

Following Wang 25 , the Ising-interaction can be treated 
in mean field approximation (MFA) by introducing a 
nearest neighbor "covalent bonding" of the Jordan- 
Wigner fermions \ = ( c l c i+i)- 

H M f = J - 2 x) cos k c\ck with x = ~- (3) 

k 

The ground state of the Heisenberg model has no net 
magnetization (Sf) = 0. Within the fermionic repre- 
sentation this implies that the fermion system is at half 
filling (cjcj) = 5. The ground state is obtained by fill- 
ing up all negative energy states. This leads to a Fermi 
surface at wave vectors kf = ±f , as displayed in Fig. 1. 
Adding/removing a fermion to the system corresponds to 
S z = ±l-excitations, S z = 0-excitations can be realized 
by particle-hole excitations. The particle-hole continuum 
of the Jordan- Wigner fermions, displayed in Fig. 1, is 
very similar to the two-spinon continuum 26 . The upper 
cutoff of the Jordan- Wigner particle-hole continuum is at 
(2 + -) J ss 3.27J and therefore close to ttJ which is the 
maximum energy for two spinons. 

B. RPA for optical conductivity — d=l 

High-energy spin excitations can be observed in the 
mid-infrared range of the optical conductivity u(lo) due 
to the simultaneous excitation of a phonon. 22 The optical 
conductivity of the Id spin-chain compound Sr 2 Cu03 24 
has been nearly perfectly reproduced by Lorenzana and 
Eder 27 using an ansatz based on numerical results in fi- 
nite chains, sum rules and Bethe ansatz results. Origi- 
nally, a similar procedure was suggested by Miiller et al. 28 
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FIG. 1. Left panel: mean field dispersion for the Jor- 
dan- Wigner fermions in d=\, grey shading denotes the fill- 
ing in the ground state. Right panel: continuum of Jor- 
dan- Wigner particle-hole excitations in d=l. 
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for the evaluation of the dynamic structure factor S(k, u), 
taking advantage of the observation that the two-spinon 
contribution is the class of Bethe-ansatz solutions which 
carries most of the weight of the continuum excitations. 
Only recently, it has been possible to determine the two- 
spinon contribution to S(k,u>) exactly. 29 ' 30 

For the optical conductivity, however, an exact expres- 
sion of the two-spinon contribution is not yet available. 
Nevertheless, the evaluation of Lorenzana and Eder, 27 
which nearly perfectly reproduces the shape of the cusp- 
like, wide structure in <j(oj), confirms that the observed 
resonance indeed results from two-spinon excitations of 
the nearest-neighbor Heisenberg model. This motivated 
us to use the established cr(u>) of the Id spin-chain as 
a reference and check for the quality of the results of 
our analytical Jordan- Wigner approach. We calculate 
the two-particle correlation function a(w) within an ex- 
tended RPA-scheme, i.e., by summing up bubble- and 
ladder-diagrams and compare our result with the exper- 
imental optical conductivity of Sr2Cu03. 24 

For the one-dimensional spin-chain the phonon- 
assisted magnetic contribution to the optical conductiv- 
ity is given by 22,27 

a(u) ~ -16^^sin 4 (|)Im(( ( 5 J B_ p ; ( 5i? p )) (L ,_ L , p , i) . (4) 
p 

The spin-flip operator SB p is expressed in MFA by 
6B p=jfH eipri (SiS i+1 - (SiSi+i)) 

i 

M e-*/*-L £(i - 2 X ) cos(fc + |) c t Cfc+p . (5) 

k 

This yields for the dynamic spin-flip correlation function 
in Zubarev notation 

((SB_ p ;SB p )) = 1 £(1 - 2 X ) 2 {cos 2 P -B^{ Pl u) 

p 

- C os P -sm P -(B^{p^) + B^{p^)) 
+ sin 2 |i?( 2 < 2 )(^)} (6) 

with particle-hole propagators 

B^")(p,u) = £/OT<4c fe+P ;4 +p c g >) (7) 

k,q 

and the following form factors 

/ fe ° = l , fl= cos k , /, 2 -sinfc. (8) 

Summing all particle-hole scattering processes, as illus- 
trated in Fig. 2 in diagrammatic terms, a simple expres- 
sion for the renormalized particle-hole propagator can be 
obtained 




FIG. 2. Diagrammatic scheme for the extended RPA treat- 
ment of Jordan- Wigner fermions in Eq. (9). 

B^ (p, w) = b ( ^ } (p, w) + 2 J cos p b ( ^ 0) (p, w)B< 0l "> (p, u) 

- 2Jb^(p,uj)B^^( P ,iu) 
-2Jb ( ^(p,uj)B (2 ^( P ,Lu), (9) 

where the nonintcracting particle-hole propagators are 
given by 

_ (1 - (n k })(n p+k } \ 
uj + e k - e p+k - i0+ J 

Evaluation of these equations determines cr(uj) which 
is shown in Fig. 3 in comparison with the experimental 
spectrum of Sr2Cu03 taken from Suzuura et al. 24 

A simple analysis of the experimental line shape of 
Sr 2 Cu03 based on Jordan- Wigner fermions has already 




Energy [eV] 

FIG. 3. o(u>) obtained with Jordan- Wigner fermions in 
comparison with the experimental optical conductivity of 
Sr2Cu03 (solid line), taken from Suzuura et al. 24 Dotted line: 
mean-field approximation, dashed line: RPA-approximation. 
Following Lorenzana and Eder 27 we have subtracted the same 
linear background from the experimental data and we have 
used the same value for the exchange coupling J=0.246eV 
and for the phonon- frequency oj p h=0.08eV. 
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been discussed by Suzuura et al. 2A in combination with 
the experimental results. However, they restricted the 
evaluation to the XY-model which corresponds to our 
mean field evaluation apart from a renormalization of the 
energy scale by a factor of 1+2 /tt in Eq. (3). We find that 
it is important to treat the two-particle correlation func- 
tion a(uj) at least within RPA. The resonance is shifted 
to lower energies compared to the mean field approxima- 
tion. In addition we observe a cusp at u> — J(l+2/ir) as a 
precursor of the logarithmic singularity found by Loren- 
zana and Eder. 27 Despite the good agreement with re- 
spect to the cusp and to the high energy side of the cusp, 
however, our RPA treatment seems to overestimate the 
effect of the interaction because too much spectral weight 
is shifted to energies below the cusp. This effect could 
possibly be compensated by consideration of higher order 
correlations. 



III. JORDAN- WIGNER TRANSFORMATION 
FOR THE TWO-LEG S = i-LADDER 

Motivated by the results of our Jordan- Wigner fcrmion 
treatment for the Heisenberg S — i-chain we "slightly 
increase" the dimensionality and extend the approach to 



* ) » ) * ) t ) i ) * 



the nearest neighbor Heisenberg two-leg S 



-ladder: 



i,iSj+i,i+ Sj^Sj+i^) , (11) 



where Jj_ is the exchange coupling along the rungs, J the 
coupling along the legs, i refers to the site index along the 
legs and the subscripts 1, 2 label the two different legs. 

Generalizations of the Jordan- Wigner transformation 
to higher dimensions have been suggested 31,32 and may 
be adopted for spin ladders. For spin ladders the one- 
dimensional Jordan- Wigner transformation can be ap- 
plied directly, when all spins are arranged in a one- 
dimensional sequence. With this scheme we can control 
the range of the interaction terms through a convenient 
choice of a path covering all sites, whereas the appli- 
cation of a two-dimensional representation to the spin 
ladders would generate long-range interaction terms in 
the Hamiltonian. Possible path configurations through a 
two-leg ladder arc shown in Fig. 4. 

The path displayed in Fig. 4(a) is obviously very close 
to the one-dimensional situation. As a consequence the 
rung interaction is difficult to treat in this representation 
because every product of neighboring rung-spins contains 
an infinite number of phase-factors. The rung-coupling, 
however, is a relevant perturbation since the excitation 
spectrum of a two-leg ladder remains gapped for all cou- 
pling ratios Jj_/ J. Therefore a path which passes through 
all the rungs should be more suitable. Possible realiza- 
tions are a zigzag path 21 and a meander path 20 , dis- 
played in Fig. 4(b) and (c), respectively. Although the 
zigzag path appears simpler and more symmetric at first 
sight, we show in Appendix A that within a mean field 
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(b) 
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FIG. 4. Possible path configurations for a two-leg ladder. 

treatment (analogous to the previous section) only the 
meander path can reproduce the correct strong coupling 
limit of the one-triplet dispersion efe = Jj_ + J cos k for 
J±/J 3> 1. Therefore we chose the meander path even 
though the nearest neighbor Hciscnberg-Hamiltonian is 
slightly more complicated in this representation. 

Following Dai and Su 20 we divide the ladder into two 
sublattices as indicated in Fig. 5. Introducing two species 
of spinless fermions a, and fa, the spin operators on the 
two sublattices transform as: 



(12) 



where the summation in the phase factor is along the 
meander path. For products of spin operators, which 
are not successive along the meander path, like e.g. 
^ta^t+i.p + S^ a S+ +l f3 , the phases corresponding to in- 
termediate sites along the meander path do not can- 
cel. This is different from the one-dimensional situation 
where all nearest-neighbor spin operators are also suc- 
cessive along the path. Using transformation (12) the 
Heisenberg Hamiltonian of Eq. (11) becomes: 

H = J±^2 {^(atft + Pj ai ) + (atai - \xj3\pi - \) J 



+ 3 E \ \ + al/W^"^"^ + H.c. 



(13) 



Unfortunately the phase factor e l7T ( n i 3 i +n °'i+i') from spin 
products of non-successive sites cannot be treated ex- 
actly. Dai and Su have replaced this phase factor by its 
average value. This treatment, however, can be improved 
in a systematic way by expanding the phase factor: 

e iAn, i+ n ai+1 ) = (1 _ 2f3 t mi _ 2a \ +1 a i+1 ) (14) 
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Si-i.B Sj a S i+1 13 

FIG. 5. Sublattice structure for the meander path. 

and reinserting this expansion into Hamiltonian (13). In 
this way we obtain additional interaction terms contain- 
ing 4- and 6-fermion operators which we now treat on 
the same footing as the Ising-interaction terms. 

A. Mean field treatment — Ladder 

Inspired by the results of the mean-field treatment for 
the spin chain, we adopt a similar approach here and 
consider all possible nearest-neighbor bond amplitudes: 



Xo 



, Xi = (Pt<*i+i) , X2 = • (15) 



Taking into account all possible contractions of the 4- 
and 6-fermion operator terms we arrive at the following 
mean-field Hamiltonian: 



H MF = ^(ikalPk + H.c.) 



(16) 



with 



lk = J±{^ -Xo) +4JxoXi (17) 
+ J cos fc(i + 2x1 - 4 XiX2 - Xi - X2 - 2x?) 



U sinfc(- + 2\l ~ 4 XiX2 



Xi +X2 + 2 X i). 



This expression has already been simplified using that xo, 
Xi and xi turn out to be real. The above Hamiltonian 
can easily be diagonalized 

Hmf = tkiaj&k - Pjfik) with e k = \y k \ (18) 



using 



-y=u k {a k + P k ) 



u k = v k 



i<p k /2 



Pk = -j=v k (a k - f3 k ) 

lk = \lk\e 1 ^. (19) 



For the isotropic ladder, i.e. J = Jj_, we obtain xo = 
-0.3617, xi = -0.2679 and X 2 = 0.1777. As the product 
of the bond amplitudes around a plaquette is negative, 
our mean field treatment corresponds to a 7r-flux state 
of the spinless fermions. Note, that this 7r-flux phase is 
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— JW - mean field 

- - JW - Dai, Su 
• DMRG 
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0.5 1 

k/n 

FIG. 6. Dispersion for the isotropic ladder J = J±. 
Solid line: mean field dispersion for Jordan- Wigner fermions; 
dashed line: dispersion obtained by averaging the phase fac- 
tor analogous to the treatment by Dai and Su; 20 circles: 
one-triplet dispersion obtained with DMRG for a A?=80-site 
ladder. 17 



different from that of Azzouz et al., 21 who have by con- 
struction replaced the complete phase factor by a 7r-flux 
phase. Both approaches are compared in Appendix C. 

The resulting dispersion is displayed in Fig. 6, in com- 
parison with the dispersion for a iV=80-site ladder ob- 
tained by DX1RG. 17 For momenta between k w 0.57T- 
0.97T we find nearly perfect agreement with the DMRG 
results. The spin gap at k — w, however, is too small 
and for momenta k < ir/2 the energy is overestimated. 
Nevertheless, our mean-field treatment is able to repro- 
duce a dip for small momenta, which is a precursor of the 
symmetric (with respect to k = ir/2) spinon dispersion 
of the spin chain. The dip becomes more pronounced 
when the leg coupling J/J± is increased as can be seen 



H 2.0 




FIG. 7. Mean-field dispersion for Jordan- Wigner fermions 
for different coupling ratios J/J±. Inset: spin gap A s as a 
function of J/J±. 
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in Fig. 7. Note that within a mean-field treatment of 
the bosonic bond-operator representation of elementary 
rung-triplets 10 ' 11 it has not been possible to achieve a 
dispersion-dip for the isotropic ladder. 

To demonstrate the improvement of our mean-field 
evaluation with respect to the mean-field treatment by 
Dai and Su, 20 who have replaced the phase factor by 
its expectation value, we have added their dispersion in 
Fig. 6. Qualitatively it is very similar to our mean-field 
dispersion. Its magnitude, however, is by a factor of 
about 1.5 too small over a large part of the Brillouin 
zone. Therefore we conclude that an adequate treatment 
of the phase factor is very important and it is necessary 
to go beyond averaging of the phase factor. 

The spin gap A s , which is shown in the inset of Fig. 7 
as a function of J/ J± , is underestimated by our approach. 
It even vanishes for J ss 1.7 J ±. As the spin gap is known 
to be finite for all finite J/Jj_, this coupling ratio marks 
the breakdown of our mean field treatment. 



B. Extended RPA-treatment — Ladder 

The optical conductivity a(u>) constitutes a powerful 
probe of the spin excitations of a spin ladder. It has 
been possible to verify experimentally the existence of 
a 5=0-two triplet bound state in the optical conduc- 
tivity of the spin ladder compound (La,Ca)i4Cu2404i 16 
and a more refined analysis has shown that it is neces- 
sary to include a 4-spin cyclic exchange interaction 17 of 
about J cyc w 0.20 — 0.27Jj_. Bound states in the singlet 
and triplet excitation channel were predicted by Sushkov 
and Kotov 11 and by Damle and Sachdev, 33 similar to 
the states analyzed by Uhrig and Schulz 34 in dimerized 
spin chains. More extensive perturbative investigations, 
within a linked cluster expansion, were performed by 
Trebst et al. 12 and the observability of the singlet bound 
state was suggested by Jurecka and Brenig 35 . 

Here, we will demonstrate, how the spin-flip corre- 
lation functions, which contribute to cr(u>), can be ob- 
tained within the Jordan- Wigner approach. We discuss 
the resultant correlation functions for an isotropic ladder 
J = Jj_, without cyclic exchange (J cyc = 0), and focus on 
the S*=0-bound state and the continuum excitations. For 
the calculation of the optical conductivity we will concen- 
trate on an isolated Cu203-laddcr. The phonon-assisted 
magnetic contribution to a(ui) results largely from the si- 
multaneous excitation of two neighboring spin-flips plus 
a Cu-0 bond-stretching phonon: 

a(uj)~-uY^ E /pM(«B-p;«Bp))(^) (20) 

P Py=0,TT 

where p = (p,p y ) and the operators 

sb p s = ^EE eiprt, ( s *.' s *+M-<Si,is i+M )) 

i 1=1,2 



«?p ung = Jf E e<PP ' ( S *.i S *.2 - (Si,lSi,2» (21) 

i 

are the spin-flip operators for polarization of the electrical 
field along the legs and the rungs, respectively. 

Following our previous treatment in Ref. 17, we con- 
sider phonon form factors given by: 

/^ = 8sin 4 (|) , /™ ng = 8sin 2 (|) + 4. (22 ) 

Here, f p eg originates from the coupling to in-phase and 
out-of-phase stretching modes of O-ions on the legs and 
is the same form factor as for an isolated spin chain. For 
/™ ng we take in addition to the out-of-phase stretching 
mode also the vibration of the O-ion on the rung into ac- 
count, which is responsible for the constant contribution 
in Eq. (22). 

For the calculation of the spin-flip correlation function 
we apply the Jordan- Wigner transformation (12) to the 
spin-flip operators S^iSj+i,! ± Sj^Sj+i^ and Sj.iS^. 
All terms with 4 and 6 fermion operators are reduced 
to two-operator terms by replacing all surplus operators 
with their contractions (15). With this procedure the 
Fourier transform of the spin-flip operators becomes: 

SB™* = J= £(1 - X o) (alP p+k + P\a p+k ) (23) 

V k 

5B>;i =0 = J= E {4/W« + be ik ce-W)) 

V k 

+ Pla p+k {a + be-'b+Q - ce*)} 
= -!= E {a\(3 P+k+ .{-a + be ik ce-^ +fe )) 

+ pla p+k+7T (-a - be-^+V + ce lk )) 

with 

a = 4xoXi , 6 = i + 2xo - 4xiX 2 - Xi 

c = 2x? + X2 , b = i - 2xo + 4xiX2 - Xi • (24) 

Inserting the spin-flip operators (23) into the optical 
conductivity (20) produces a sum of particle-hole prop- 
agators with different form factors. To evaluate these 
particle-hole propagators in RPA we prefer to use the 
original fermionic operators a, (3 because transformation 
to the operators a, (3 (19) would increase the number of 
interaction terms considerably. 

Prior to the derivation of the RPA-equations, the in- 
teraction terms in the Hamiltonian have to be reduced 
to two-particle interactions in order to deal only with 
4-particle vertices. Accordingly, all 6-operator terms, 
which appear in Eq. (13), are reduced to 4-operator 
terms by replacing all possible contractions with the cor- 
responding bond amplitudes (15). In this way, we ob- 
tain the following reduced interaction term from Hamil- 
tonian (13): 
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FIG. 8. Diagrammatic representation of the processes 
which contribute to the extended RPA-treatment, in real 
space. 



±L + e Kk 2 -k 3 ) + + 4 X2 ) e -«(fe2-fc3) 



- 2xo(e 2fel + e- lk2 + e lks + e^ 4 )] 

+ 4,4,^3^4 2 X ie 4 ( fel - fe4 )} (25) 

A set of RPA-equations for the particle-hole propaga- 
tors, which are listed in Appendix B, can be obtained 
by consideration of all possible vertex-configurations of 
the interaction Hamiltonian (25). In real space these 
vertices correspond not only to a summation of bubble- 
diagrams, but also include ladder-diagrams and other 
non-local terms as indicated in Fig. 8. For this reason 
we use the term extended RPA-treatment. 



C. Spin-flip correlation functions 

We have calculated the correlation functions for spin- 
flips along the legs SB l °^ =0 , 6B l p^ =7T and for spin-flips 
along the rungs SB™" 6 for an isotropic ladder J = J± 
using the RPA-treatment of the previous section. The 
results are presented in Figs. 9, 11 and 12, respectively. 
For comparison the mean-field evaluation of each of the 
correlation functions is displayed in the lower panels. 

In the mean- field evaluation of ((SB lc ^„ _ n ;SB 1 ^ _ n )) 
(lower panel of Fig. 9) one observes van-Hove singular- 
ities at the upper edge of the continuum for small mo- 
menta and at the lower edge of the continuum for large 
momenta. With the RPA-treatment the van-Hove sin- 
gularities at the continuum edges disappear. For small 
momenta the maximum of the continuum is shifted from 
the upper edge downwards to about u> sa 3 J. At large 
momenta we see the formation of the S'=0-bound state. 
The bound state emerges from the continuum at ap- 
proximately k 0.37T, it passes through a maximum at 
k ss tt/2 and a minimum at k = n. In Fig. 10 the dis- 
persion of the bound state is compared with the DMRG 
calculation for a A^=80-site ladder. 17 We find good agree- 
ment between both methods, only the energy of the RPA- 
dispersion is slightly too low. This indicates that the 
interaction strength is somewhat overestimated in RPA. 



In the mean field evaluation of the out-of-phase com- 
ponent of the correlation function for spin-flips along the 
legs ((SB^^iSB**^)) (lower panel in Fig. 11) the 
van-Hove singularities at the continuum edges are sup- 
pressed. The overall momentum dependence on p ap- 
pears to be reversed when the out-of-phase component in 
Fig. 11 (upper panel) is compared to the in-phase com- 
ponent in Fig. 9 (upper panel). This is caused by the 
checker-board sublattice structure of the meander path, 
which shifts the momentum of the particle-hole propa- 
gator 6B^& _ n by 7T in relation (23). The inversion of 
the momentum dependence is especially noticeable for 
the bound state. However, the out-of-phase component 
should not contain the bound state at all but only con- 
tribute to the continuum excitations, an issue that we 
have addressed previously in Ref. 17. The argument is 
based on the observation that the out-of-phase compo- 
nent originates from the excitation of 3 different rung- 
triplets, 1 when it is expressed in terms of rung-triplet 




co /J 

FIG. 9. RPA and mean-field evaluation of the momentum 
resolved correlation function ((5B_fS p =0 ;5B^ _ )), where 
a broadening of S = 0.01J has been used. The grey lines 
(upper panel) designate the S=0-bound state. The dashed 
dark grey line is a projection of the dark grey line in the lower 
panel, which links the points of sharp increase in the middle 
of the continuum. This is a precursor of the upper edge of 
the 2-spinon continuum in single chains. 23,36 In RPA a dip 
structure remains at the same position in the continuum. 
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FIG. 10. Particle-hole continuum of the Jordan-Wigner 
fermions and S'=0-bound state (thick solid line) in com- 
parison with the two-triplet continuum (open symbols) and 
5=0-bound state (filled symbols) obtained by DMRG for a 
ladder with 7V=80-sites. 17 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



RPA 



s u S i+i,r S i,2 S i+i,2 




Mean Field 




FIG. 11. RPA and Mean-field evaluation of the correla- 
tion function {(SB 1 ^ =7r ; SB 1 p t f y=n )) , where a broadening 
of S — 0.01J has been used. The grey lines indicate the 
S^O-bound state. 




FIG. 12. RPA and Mean-field evaluation of the correlation 
function ((SB T ™ S ; SB r p uns )), where a broadening of S = 0.01J 
has been used. The grey lines indicate the 5^0-bound state. 
The dashed dark grey line in the upper panel is the same as 
in Fig. 9. 



operators. 9 ' 10 The S'=0-bound state, on the other hand, 
arises from scattering processes of two equal triplets and 
therefore cannot be present in the out-of-phase compo- 
nent. Although the spectral weight of the bound state 
in the p y — 7r component is only small, this indicates a 
failure of our approach. This failure is probably related 
to the fact that the meander path breaks certain sym- 
metries of the original ladder model. In our mean-field 
evaluation we find, e.g., (S^iSi+1,1) ^ (S^Si+i^)- This 
artificial symmetry breaking will be discussed in more de- 
tail in Appendix D. However, the resultant effects may 
average out to some extent in the sum of the spin-flip 
operators on both legs S^iSi+^i + S^. 28^+1,2, whereas 
they will probably even be enhanced in the difference 
Si.iSi+1.1 — Si^Si+i.2- This may explain the shortcom- 
ings of our approach with respect to the out-of-phase 
component, while we find reasonable results for the in- 
phase-component. 

In Fig. 12 the RPA- and mean-field evaluation for 
the correlation function of spin-flips along the rungs 
((S B 1 ^ 1 ^ ; S Bp Un& )) are shown in the upper and lower 
panel, respectively. For p = it is identical to the 
p = 0-component of the in-phase correlation function for 
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spin-flips along the legs {(SB lcs 

correspond to the correlation function for the Raman re- 



r Py ^ SB p!py=0)) aild b ° th 



37.38 



sponse."-" For larger momenta, the spectral weight of 
the rung-correlation is much smaller than the in-phase 
component of the leg-correlations and at p = ir the weight 
of the S=0 bound state vanishes according to a selection 
rule. 16 



D. Optical conductivity 

Once the momentum dependent spin-flip correlation 
functions, shown in Figs. 9, 11 and 12, are known the op- 
tical conductivity <t(oj) in Eq. (20) can easily be obtained 
by integration. In Fig. 13 the momentum integrated 
spin-flip correlation functions are displayed whereby we 
present them with form factors equal to unity (dashed 
lines) and with the form factors of Eq. (22) (solid lines) . 
To demonstrate their contribution to a(uj), all spectra 
have already been multiplied with frequency w. 

For polarization parallel to the legs and p y = (top 
panel of Fig. 13) two dominant peaks appear at w\ ~ 1.4 J 
and Ld 2 ~ 2.0J. They are caused by van Hove singular- 
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FIG. 13. Momentum integrated spin-flip correlation func- 
tions, which contribute to cr(w) for polarization along the legs 
(top and middle panel) and for polarization along the rungs 
(bottom panel). In order to visualize the resulting contribu- 
tion to a{u>) in Eq. (20) each correlation function has been 
multiplied with the frequency ui and for each of them the 
momentum integration without prefactor (dashed lines) and 
with the corresponding prefactors of Eq. (22) (solid lines) are 
shown. 
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FIG. 14. Comparison of the Jordan-Wigner spectra for the 
isotropic ladder using an exchange coupling of J— 1 100cm - 1 
with the experimental u{u>) of La5.2Ca8.8Cu24O.41. The 
experimental spectra have been measured by Windt and 
Gruninger. 17 For the phonon a frequency of cj p h=600cm -1 has 
been assumed for polarization along the legs and w p h=650 -1 
for polarization along the rungs. 



ities arising from the dispersion of 5=0-bound state 16 
at p = 7T and p w tt/2. The upper peak at u>2 is 
suppressed by consideration of the relevant form factor 
sin 4 (p/2). Also the continuum excitations are reduced 
considerably because they originate largely from small 
momenta, where no bound state is present, see Fig. 9. 
The out-of-phase component (middle panel of Fig. 13) is 
an order of magnitude smaller than the in-phase compo- 
nent. We consider the results for the out-of-phase compo- 
nent, however, less reliable due to symmetry breaking by 
the meander path, as mentioned in the previous section. 

For polarization parallel to the rungs (bottom panel of 
Fig. 13) two form factors contribute to <j(oj), a constant 
and a sin 2 (p/2)-term. Only the upper bound state at u>2 
is present in the momentum integrated spectrum because 
the S'=0-bound state is suppressed at p = n in accord 
with a selection rule, 16 see also Fig. 12. 

In Fig. 14 our spectra are compared with the exper- 
imental optical conductivity of La 5 .2Ca 8 . 8 Cu2404i for 
polarization along the legs (top panel) and for polar- 
ization along the rungs (bottom panel). Following our 
previous treatment in Ref. 16, an exchange coupling of 
J=1100cm -1 , a phonon frequency of cjp/j^OOcm -1 for 
polarization along the legs and Wj,/j=650cm -1 for polar- 
ization along the rungs have been used. The difference of 
50cm _1 in the phonon frequencies for polarization along 
the legs and the rungs accounts for the shift observed 
experimentally with respect to the upper bound state in 
both polarizations. 16,17 

For polarization along the legs we find good agree- 
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ment with respect to the bound states. The continuum 
contribution, however, is strongly underestimated. To 
some extent this is due to the fact that, so far, we con- 
sidered only the creation of two fermions at the exter- 
nal current vertex: we approximated the spin-flip op- 
erators in Eq. (23) by a particle-hole creation opera- 
tor. Without this approximation the correlation func- 
tions ((SB l ^ py=Q y,5Bp^ y=0 ^)) would contain also the 
excitation of 4- and 6-fermions, which would increase 
the amount of high-energy excitations. The weight of 
these higher order excitations will be estimated in the 
next section. 

For polarization parallel to the rungs we find very 
good agreement even with respect to the weight and line 
shape of the continuum. Rung correlations are obviously 
treated quite accurately in our meander-path represen- 
tation. There are no phase factors in the product of 
two spin operators on the same rung because they are 
neighboring along the meander path. Therefore the only 
4-fermion process results from the Ising-like term in the 
rung operator (21). 



E. Higher order contributions 

The contribution of 4- and 6-fermions to the spin-flip 
correlation functions in leg polarization can be obtained 
analogous to the 2-fermion contribution. In the spin-flip 
operator <5i?p 0g , Eq. (21), the spin operators are replaced 
by fermionic operators. Terms with 6 operators are des- 
ignated as the 6-fermion part. Terms with 4 and 6 oper- 
ators contribute to the 4-fermion part, when all surplus 
operators are replaced by their contractions, analogous 
to Eq. (23). With this scheme one obtains for the 4- and 
6-fermion part of the spin-flip operator 
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for the in-phase component of the leg-polarization. The 
out-of-phase component is very similar, only the mo- 
mentum p has to be replaced with (p + 7r), and in 
the term \ (e 1 ^ 2 "^) + e^- k ^)a\^\^ k3 a ki of the 4- 




FIG. 15. Mean-field weights for the spin-flip correlation. 
They correspond to the evaluation of one, two and three non- 
interacting particle-hole lines for the excitation of 2-, 4- and 
6-fermions displayed in (a), (b) and (c) respectively. 



fermion component the parentheses have to be replaced 
with ( e i (' C2 ~' C3 ) — e i ( k i- k i)y 

In order to estimate the contribution of processes in- 
volving the excitation of a different number of fermions, 
we compare the weights 



K y = - 



p 



■> SB p,pJ)- 



(27) 



which are obtained by integrating the n=2, 4 and 6 
fermion part to the leg-correlation functions over fre- 
quency u> and momentum p using form factors f p = 1 and 
f p = sin 4 (p/2). To keep this evaluation as simple as pos- 
sible the correlation functions are evaluated in mean-field 
theory, i.e., they are replaced by noninteracting particle- 
hole lines as indicated in Fig. 15. The frequency integrals 
in Eq. (27) can be eliminated using: 



/>oo -i />oo 

n n a i = ~ im / duj eiuj0+ n o - / dn i x 

j JO j J_ oc 



(28) 



Oj - < j + /'I I % + ej - '(I 



where a,j/(Qj =p €j ± i0 + ) are general expressions for re- 
tarded and advanced Greens functions. The remaining 
momentum integrals can be easily evaluated numerically. 

The resulting weights of the 2-, 4- and 6-fermion pro- 
cesses are displayed in Table I. The 4- and 6-fermion 
processes make up only 24% (15%) of the in-phase part 
of the spin-flip correlation function when a form factor 
of f p = 1 (f p = sin 4 (p/2)) is used. On the other hand, 

TABLE I. Mean-field weights of the 2-, 4- and 6-fermion 
contributions to the in-phase (p y = 0) and the out-of-phase 
(p y = 7r) component assuming noninteracting particle-hole 
lines. For the 2-fermion contribution the weight of the 
RPA-evaluation is given in parentheses. 
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they generate a major contribution to the out-of-phase 
component, i.e., 65% (50%). Note, that these values re- 
fer only to noninteracting particle-hole propagators and 
may be changed when the evaluation is improved. For 
example the weight of the 2-fermion contribution to the 
in-phase-component is enhanced when the particle-hole 
propagator is evaluated in RPA. In Table I these RPA- 
weights are added in parentheses. 

Since the 4- and 6-fermion processes contribute only 
to the continuum excitations, they will increase the high- 
energy weight and therefore improve the consistency with 
the experimental spectrum for polarization along the legs 
in Fig. 14. However, the consideration of 4- and 6— 
fermion processes will presumably not suffice to account 
fully for the large high-energy weight observed experi- 
mentally. One reason is, that the investigated spin-ladder 
compound La5.2Ca8.sCu2404i corresponds not exactly to 
an isotropic spin ladder J = J± but should rather be 
modeled by J 1.3J± including a finite 4-spin cyclic ex- 
change interaction of J cyc w 0.2Jj_. 17 The presence of a 
finite cyclic spin exchange J cyc and the anisotropy of the 
exchange couplings J/Jj_ > 1 will increase the amount of 
high-energy continuum excitations. 39 On the other hand 
it might also be necessary to include higher order ver- 
tex corrections. They may shift some weight from the 
bound state to the continuum excitations, because the 
RPA overestimates the interaction strength as discussed 
in Sec. IIIC. They may also help to restore the symme- 
tries which are broken artificially by the meander path 
as discussed in Appendix D. 

IV. CONCLUSIONS 

In this paper we have discussed an approach to obtain 
dynamic correlation functions in low dimensional quan- 
tum spin systems based on the Jordan- Wigncr transfor- 
mation. We have shown, that the application of standard 
perturbation theory to the new fermionic operators is rea- 
sonable even in the high-energy range. We have tested 
our approach by calculating the dynamical spin-flip cor- 
relation function for the Id spin chain, which corresponds 
to the phonon assisted magnetic contribution to the opti- 
cal conductivity, and we have found good agreement with 
the experimental spectrum of Sr2Cu03. 

Using a meander-path like arrangement of the spin op- 
erators we have extended this approach to the two-leg 
S = ^-ladders. In contrast to the Id spin chain, how- 
ever, phase factors remain in the Hamiltonian. Expand- 
ing these phase factors creates new interaction terms, 
which can be treated within standard perturbation the- 
ory. We have calculated the one-particle dispersion in 
a mean-field treatment based on nearest-neighbor bond 
amplitudes and the optical conductivity within an ex- 
tended RPA-scheme. 

For polarization along the rungs we find good agree- 
ment with the experimental spectrum of the spin-ladder 



compound La5.2Ca8.sCu2404i. Spin-flip processes on the 
rungs are represented quite reliably in our approach, be- 
cause two spins on the same rung are also neighboring 
along the meander path and therefore no phase factor 
is present in the corresponding spin-flip operator. For 
spin-flip processes along the legs, however, a phase fac- 
tor appears, which results in the excitation of 4- and 
6-fermions in addition to the 2-fermion processes. Eval- 
uation of the 2-fermion contribution to the spin-flip cor- 
relation function for polarization along the legs yields 
good results with respect to the S^O-bound state. The 
weight of the high-energy continuum, however, is under- 
estimated. To some extent this can be compensated by 
the consideration of 4- and 6-fermion excitations. 

In principle our approach can be generalized straight- 
forwardly to spin-ladders with more than two legs. How- 
ever one should think about a less refined treatment of 
the phase factors in order to reduce the increasing num- 
ber of interaction terms resulting from the longer range 
of the phase factors. 
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APPENDIX A: STRONG COUPLING LIMIT 

In this appendix we demonstrate that our mean-field 
treatment for the Jordan- Wigner fermions, which is 
based on the meander path, can reproduce the correct 
strong coupling limit = J± + Jcosfc. This is not true 
for the zigzag path within a straightforward extension of 
our mean-field approach. 

1. Meander-path 

Expanding around the rung-dimer limit we obtain: 

• zeroth order: (J/Jj_)° 
For J = the off-diagonal part of the mean-field 
Hamiltonian (16) reduces to 

lk = J±(l-Xo). (Al) 



The resulting Hamiltonian can be diagonalized eas- 
ily using Uk — Vk — 1 in Eq. (19) which yields for 
the nearest neighbor bond amplitudes 




k 
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This is the limit of rung dimers and the correct 
value for the energy of a single rung-triplet excita- 
tion is obtained as 



with 



£fe = hk\ = J± ■ 



(A3) 



• first order: (J/Jl) 1 
Re-substituting the zeroth order bond amplitudes 
in Eq. (17) yields 7^ = J± + J(cos k + i sinfc). Ex- 
panding the resulting dispersion to first order in 
J I Jj_ one obtains the correct strong coupling limit: 



£fe = \lk\ = J± 

~ J± + J cos k . 



l + 2^-cosfc+(^-) 2 
J± J± 



(A4) 



This expression contains already all terms to or- 
der J I J_l . This can be easily seen when the above 
results are inserted into Eq. 15 for the bond am- 
plitudes. To order J/J± one obtains for the di- 
agonalization transformation Uk — Tfe / 1 Tfe | ~ 1 + 
iJ/J± sin k and consequently for the bond ampli- 
tudes xo = -1/2 and xi = ~X2 = J/(4J±). 
Whereas \o remains unchanged, Xi and X2 are pro- 
portional to J/Jj_ and therefore contribute to 7^ 
and €k only in second order. 



2. Zigzag— path 

In this subsection we discuss the mean field treatment 
for the zigzag path. In particular we show that it is does 
not reproduce the correct strong coupling limit. 

For an illustration of the zigzag path see Fig. 4(b). 
Here, the underlying sublattice structure is very simple, 
each leg corresponds to a sublattice. Choosing, e.g., in 
the Jordan- Wigner transformation, Eq. (12), 5j iQ on the 
lower leg and Si^ on the upper leg of Fig. 4(b), the XY- 
part of Hamiltonian along the legs becomes 



H 



XY.leg 



a\a i+1 (l-2p}p i ) (A5) 
+ /3t/3 i+1 (l-2at +1 a i+1 )+H.c.]} . 



Here the expansion of the phase factor yields only 4- 
operator terms. Taking into account all nearest and next 
nearest neighbor bond amplitudes: 



Xo = {0-ai) , xi = (ajai+i) 
X2 = (Pjfii+i) , Xs = (Pjoti+i) , X4 = (atfii+i) 

one obtains the following mean-field Hamiltonian 

Hmf = ^2 \jo,ka\Pk - Ja,kO:{ak ~ H3,k0[Pk + H.c.} 
k 

(A7) 



(A6) 



lo.k = J±{^ - Xo) + 2J( X ae lk + X s) 
la.k = Jxie lk , l(3.k = JX2e tk . 



(A8) 



For J = this equals the mean-field Hamiltonian for the 
meander path, Eq. (16), and therefore the same values 
for the bond amplitudes xo = — 5 , Xi = X2 = X3 = are 
obtained. Reinserting these values into Eq. (A7) yields 

7o,fc = J_l - J (cos k + i sin k) , -y a . k = jp ik = . (A9) 

The corresponding dispersion is to first order in J/ J± 



(,, - -, u . = 1 2-j^-cosfc • -J-' A 

J± — J cos k 



(A10) 



which obviously differs from the correct strong coupling 
limit (A4). 



APPENDIX B: RPA EQUATIONS 

The RPA-equations for the particle- hole propagators of 
the spinless Jordan- Wigner fermions can be obtained by 
considering all possible vertex configurations of the RPA- 
Hamiltonian (25). The explicit form for the renormalized 
particle-hole propagators is: 

2 



A=0 



+ 7 ta,«/3t D /3\p,atr* + 7 t/3,,5Qt n ^tp.ao-t 

(Bl) 

For simplicity the momentum and frequency indices p 
and co have been omitted. The noninteracting particle- 
hole propagators b and 6, including the appropriate form 
factors from the internal interaction vertices, are defined 
as: 

K'W = Jf E 6 7+P,fc+ (P' fc ' w ) 
fc 

- 77 E fk9 (S ' X) (P. *) fy^t (P, *, -) (B2) 
fc 

with form factors 

/£ = 1 , /^e* , fi = e~ ik , (B3) 

and 
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ff d,0) 


(P, fc ) 


= -J ± .f k +2J XQ (l + e-*n.f 2 k 




,(M) 


(P, fc ) 


= J [2 X0 (1 + e 4p )A° - axie^A 1 - (1 + 4 X2 )/ fe 2 ] 




(P, k) 


= — J (A 1 + 2 Xl e-^.f fe 2 ) 




g (2.0) 


(p, k) 


= [J ± + j(e-* + (l + 4 X 2)e <p )] A° 
-2J X oe ip fl-2J Xo fi 




,(2,1) 


(p, k) 


= 2J(-xoe ip f fe ° + xie <p / fe 1 ) 




ff (2,2) 


(p, k) 


-2J(- X0 A° + Xie -^A 2 ) 




ff (3,0) 


(P, fc) 


= [Jx + J(e* + (l + 4 X2 )e- <p )] A° 
-2J Xo fl-2Jxoe- ip fZ 




5 (3,D 


(P, k) 


= 2J(- X0 f k+X ie* p f 1 k ) 




,,(3,2) 


(P, fc) 


= 2j(-xoe- ip f k + X ie- ip fi) 




5 (4,0) 


(Pj 


= -J ± f% + 2Jxo(l + e ip )fk 




S ™ 


(P, fc) 


= -J {2 Xl e^fl + fl) 




5 (4,2) 


(P, fc) 


= j[2 Xo (l + e -^)A°-(l + 4 X2 )/ fc 1 
-2 X ie- 4p .f fc 2 ] 




ff (5,0) 


(P, k) 


= . 9 ^°)(p,A : ) = 2J Xl (e^ + e^)/ fc 






(p, k) 


= <7 (6,1) (p,fc) = -2Jxi/ fc 2 




5 (5,2) 


(P, A;) 


= 5 ( 6 ' 2 )(p,fc) = -2J X i/ fe 1 


(B4) 



The noninteracting particle-hole propagators b° are given 
as: 



^0 



ft 



ft 



"at/3,aat 



i ( 5 ?(p,A ; ,a;)- 52 (p,fc, W )) 

6 /Jta,A8t = J U P+fc (ffl(P,*,w)+ff§(p,fc,w)) 
& at/3,/3/3t = j4+k (g°i(P,k,Uj)+g ( >(p,k,Lj)) 



h° 

°/3t/3,a/3t 



u\ (g1(p,k,u) + g%(p,k,u)) 



b° 

b % n .a„i = h %,0^ = -J v l {gi(p,k,u>)+g%(p,k,uj)) 



h° 

at/3,c,/3t 
°/3t a ,/3at 
°at/3,/3at 
°/3t Q ,Q/3t 

with 



= -J U k V l+k (ffl (P. fc > - fif° (P. fc, W)) 



VfcWp+fe (ffi (P, fc > - fc, w)) 



-j^p+k {gi(p,k,uj) - g°(p,k,u)) 
-\ u l u l+k {gi{P,k,oj) -gl{p,k,w)) 



(B5) 



g°(p,k,u>) = 
gl{p,k,u) = 



l 



w — efe — e p +fe + i<5 
1 



W + £ fc + £p+fc 



(B6) 



Here, X o, X i and X 2 are the bond amplitudes (15), u k , 
v k are the coefficients for the diagonalization (19) of the 



mean-field Hamiltonian (16) and e k is the mean-field dis- 
persion (18). 



APPENDIX C: ROLE OF THE PHASE FACTOR 

In this appendix we briefly comment on the role of the 
"phase factor" in the mean-field evaluation of the triplet 
dispersion. The phase factor e J7r (™ft+™°.+i) i n Hamilto- 
nian (13) was generated by products of spin operators 
with site labels not in sequence along the one-dimensional 
meander path. As the "matrix elements" of the phase op- 
erator are ±1 one might speculate to find a reasonable 
mean field result by replacing the operator uniformly by 
±1. This corresponds to a flux phase treatment of the 
phase factor, where the flux through a plaquette is chosen 
to be and n, respectively. Applying the same kind of 
mean field treatment as in Sec. Ill A one obtains e k = \jk\ 
for the mean field dispersion, where 



7k = J±{-^~Xo) + J cosfc(l-Xi- X2 ) + iJ sinfc( X2 - X i) 
corresponds to zero flux and 



7k = J±{ 2~Xo) - J cosk(xi+ X2 ) + U sinfc(l - X i+ X2 ) 

to a 7r-flux phase. The resulting dispersions are displayed 
in Fig. 16 together with the mean field dispersion of the 
meander-path, where (thick solid line) the correlations 
related to the phase factor are taken into account and 



JW - meander path 

JW - Dai, Su 

■ ■ ■ ■ JW - 0-flux 

- - JW - it-flux 

• DMRG 




FIG. 16. Dispersion obtained for a flux-phase approxi- 
mation of the phase factor with 0-fiux (dotted line) and 
7r-flux (dashed line) through a plaquette. Thick solid line: 
mean-field dispersion obtained for the meander path in 
Sec. Ill A where the correlations related to the phase fac- 
tor have been taken into account. Thin solid line: dispersion 
analogous to Dai and Su, 20 where the phase factor has been 
replaced by its average value. Symbols: one-triplet dispersion 
obtained by DMRG for a A=80-site ladder. 17 
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where (thin solid line) the phase factor has been replaced 
by its average value. For comparison also the DMRG 
results for a iV=80-site ladder 17 are displayed. Obvi- 
ously, the form of the dispersion improves considerably 
when the phase factor is treated in a more adequate way. 
The simple replacement by a flux phase shows only poor 
agreement with the DMRG result. It is interesting to 
note, however, that the spin gap remains finite for all 
Jj_/J > within the 7r-flux treatment of the phase fac- 
tor, as has been observed by Azzouz et al. 21 The mean 
field evaluation of the phase factor (thin solid line) im- 
proves the form of the dispersion at least qualitatively. 
For a reasonable quantitative agreement with the exact 
dispersion, however, it seems to be necessary to consider 
also the correlations related to the phase factor (thick 
solid line). 



is necessary in order to obtain a consistent description of 
the experimental data for cuprate spin ladder compounds 
like (La,Ca)i4Cu2404i. 17 The realistic value of J cyc w 
0.2 — 0.27Jj_, however, is not strong enough to induce 
a static dimerization pattern. Nevertheless the correla- 
tion length for dimerization fluctuations increases with 
cyclic spin exchange and the dip in the one-triplet disper- 
sion is reduced. 17 This could to some extent explain the 
good agreement we find for our Jordan- Wigner calcula- 
tion with the experimental spectra of (La,Ca)i4Cu24C>4i, 
as our artificial symmetry breaking seems to imitate some 
of the effects of a cyclic spin exchange. 



APPENDIX D: SYMMETRY BREAKING 

A problem of applying the Jordan- Wigner transforma- 
tion by placing a path throughout the ladder is that both 
the zigzag and the meander path break symmetries of 
the original ladder configuration. This problem has al- 
ready been mentioned in the discussion of the out-of- 
phase spin-flip correlation function in Sec. IIIC. The 
sublattice structure underlying the meander path, see 
Fig. 5, restricts the translational symmetry along the legs 
to translations of an even number of sites. Furthermore, 
it breaks the reflection symmetry with respect to the 
center line crossing all rungs perpendicularly. In MFA 
this gives rise to different values for the bond amplitudes 
Xi = (fitai+i) = -0.2679 and X 2 = (at/? i+1 ) = 0.1777. 
Even though the fermionic bond amplitudes have no di- 
rect physical meaning, these discrepancies translate into 
different expectation values for the products of neighbor- 
ing spin operators. When the spin operators are trans- 
formed to fermionic operators according to Eq. (12) and 
all contractions are replaced by their mean-field values 
as before, we obtain for the product of the z-components 
(Sf, a S? +h0 ) = -|X2| 2 , (StpSf +ha ) = -| Xl | 2 and for 
the xy-components ^{S^~ S~ +1 a + S^pSf +1 a ) = xi and 

h( S t<* S «i,f> + S ta S &if) = 4 Xi(x§ - X1X2). Therefore 
two symmetries are broken with respect to the spin op- 
erators: (i) translational symmetry by an odd number 
of sites as discussed before (Si, a Si+i t p) =/= (£>i+i )i g/Si+2, a ) 
and (ii) rotational symmetry in spin space (Sf ia S? +1 p) ^ 

The physical interpretation of breaking the trans- 
lational symmetry is quite obvious. (Si tCt Si + i,p) ^ 
{Si+i,pSi+2,a) signifies that a static dimerization pattern 
is induced along the meander path. Its amplitude, how- 
ever, is very small (Si+i,pSi+2, a - S^aSi+i^) = 0.037. 
A staggered dimerization pattern of this kind can, e.g., 
be induced by consideration of a 4-spin cyclic exchange 
interaction 40 of about J cyc /J w 1 — 2. Indeed, it has 
been found that the inclusion of a cyclic spin exchange 
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